import sympy

x=sympy.Symbol('x')
f=x**2+3*x
# 一阶导数
res=sympy.diff(f,x)
# 二阶导数
res2=sympy.diff(f,x,2)
sympy.pprint(res)
sympy.pprint(res2)
